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Abstract 

The recently proposed universal relations between the moments of the polydis- 
persity distributions of a phase-separated weakly polydisperse system are analyzed 
in detail using the numerical results obtained by solving a simple density func- 
tional theory of a polydisperse fluid. It is shown that universal properties are the 
exception rather than the rule. 
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Many, natural or man-made, systems are mixtures of similar instead of identical 
objects. For example, in a colloidal dispersion |]J the size and surface charge of the 
colloidal particles are usually distributed in an almost continuous fashion around some 
mean value. When this distribution is very narrow the system can often be assim- 
ilated to a one-component system of identical objects. Such a system is usually 
called monodisperse whereas otherwise it is termed polydisperse. Since polydispersity 
is a direct consequence of the physico-chemical production process it is an intrinsic 
property of many industrial systems. Therefore, many author H have included polydis- 
persity into the description of a given phase of such systems. More recently, a renewed 
interest can be witnessed for the study of phase transitions occuring in weakly polydis- 
perse systemsQ. The phase behavior of polydisperse systems is of course much richer 
than that of its monodisperse counterpart. It is also more difficult to study theoret- 
ically, essentially because one has to cope with an infinity of thermodynamic coexis- 
tence conditions Q. Therefore, several authors have proposed approximation schemes || 
which try to bypass this difficulty. In the present study we take to opposite point of 
view by solving numerically the infinitely many thermodynamic coexistence conditions 
for a simple model polydisperse system. On this basis we have studied the radius of 
convergence of the weak polydispersity expansion used in ref.4 and found that their 
"universal law of fractionation" and some of their conclusions have to be modified in 
several cases. 

The statistical mechanical description of a polydisperse equilibrium system is equiv- 
alent to a density functional theory [|| for a system whose number density, p(r,a), de- 
pends besides the position variable r (assuming spherical particles) also on at least 
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one polydispersity variable a (which we consider to be dimensionless) . Such a the- 
ory is completely determined once the intrinsic Helmholtz free-energy per unit vol- 
ume, f[p], has been specified as a functional of p(r, a) (for notational convenience 
the dependence on the temperature T will not be indicated explicitly). For the spa- 
tially uniform fluid phases considered here (and also implicitly in ref.4) we have, 
p(r, a) — > jo(cr), and the pressure can be written as, p[p] = J dap(a)fi(a; [p]) — f[p], 
where p(a; [p]) = 8f[p]/5p(a), is the chemical potential of "species" a. When a par- 
ent phase of density po(a) phase separates into n daughter phases of density pi(a) 
(i = 1, . . . , n) the phase coexistence conditions imply that, p[pi] = p[p2\ = ■ ■ ■ = p[pn], 
and p{a\ [pi]) = p(c; [P2]) = . . . = p(a; [p n ])- For simplicity we consider here only the 
case of two daughter phases (n = 2) and rewrite moreover pi(a) = pihi{a){i = 0, 1, 2) 
in terms of the average density pi and a polydispersity distribution hi(a) such that 
/ dahi(a) = 1. Since the ideal gas contribution to f[p] is exactly known [|| one has, 
p(a; [p]) = k bT ln{ A 3 (a) p(a)} + p e x(&] [p]), where ks is Boltzmann's constant, A(<r) is 
the thermal de Broglie wavelength of species a and p ex the excess (ex) contribution to 
p. This allows us to rewrite the equality of the chemical potentials of the two daughter 
phases as, hi(a) = h,2(cr)A(a), where A(a) is a shorthand notation for: 

P2 

= — exp (3 {p ex (a; [p 2 ]) - Pex(<r] [pi])} (1) 
Pi 

with = 1/ksT. The polydispersity distributions are further constrained by the 
relation, x\hi(a) + X2h2(c) = /io(<r), which expresses particle number conservation. 
The number concentration of phase 1, x\ = 1 — is given by the lever rule: x\ = 
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__£i — . po P2 Combining these two relations one finds: 

P\—p% po ° 

h 2 (a)-h 1 (a)=h (a)-H(a) (2) 

where H(a) = (1 — A(a))/(x2 + x\A(a)). Eq.(2) is the starting point to relate the 
difference between the moments of the daughter phases, A^ = J da c fc (/12(c) — /11(a)), 
to the moments, £fc = / daa k ho(a) (k = 0, 1,2, . . .), of the parent phase distribution 
ho(a). Indeed, when a is chosen such that ho(a) tends to the Dirac delta function 5(a) 
in the monodisperse limit, A& can be obtained from (2) by expanding H(a) around 
<T = 0, H(a) = Y^i^o a i°"'> yielding for a weakly polydisperse system, A k = Y<hLo a iCi+k- 
The normalization of the hi(a) (i = 0, 1, 2) implies Ao = 0, £0 = 1 or ao = — J2f^i a z6) 

and eliminating ao from A& yields the general moment relation: 

00 

A fc = ai£k+i + ^2ai(£k+i (3) 
1=2 

where we took moreover into account that a can always be chosen such that £1 = 0. 
When only the first term in the r.h.s. of (3) is retained we recover the universal law 
Afc/A; = £fc+i/£z+ii put forward in ref.4. The question left unanswered by the study 
of ref.4 concerns the radius of convergence of the weak polydispersity expansion (3). 
In order to study this problem in more detail we now consider a simple model system 
for which we can determine the hi(a)(i = 1,2) numerically and compare the results 
with (3). The free energy density functional chosen here corresponds to a simple van 
der Waals (vdW) model0 for the liquid-vapor transition in polydisperse systems of 
spherical particles of variable size: 

f[p]= k B Tfdap(a){H^^l)-l} 

+ydaJda'V(a,a')p(a)p(a') (4) 
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where, E[p] = 1 — J dav(a)p(a), describes the average excluded volume correction for 
particles of radius R a and volume v(a) = ^-R&, while V(a,a') = J drV(r;a,a r ) is the 
integrated attraction between two particles of species a and a', for which we took the 
usual vdW form, V(r; a, a') = —eo(R a + i? CT /) 6 /r 6 for r > R a + R a > and zero otherwise, 
eo being the amplitude of the attraction at the contact of the two particles. The size- 
polydispersity can be described in terms of the dimensionless variable, a = R a /R — 1, 
with R the mean value of R a in the parent phase, hence £1 = J da aho(a) = 0. The 
thermodynamics is given in terms of ho(a), the dimensionless temperature t = ksT/eQ 
and the dimensionless density rj = Vop, with vq = -f-i?o an d Ro the value of R a in 
the monodisperse limit. The coexistence conditions are integral equations which can 
be solved numerically using, for instance, an iterative algorithm^] for any t, r/o = i>oPo 
and ho(a). For ho(a) we took a Schulz distribution^ with zero mean. The normalized 
distribution is given, for — 1 < a < oo, by ho(a) = a a (l + a) a ~ 1 e~ a ( 1+a ^ /T(a), with 
r(a) the gamma function and 1/a a width parameter which measures the distance to 
the monodisperse limit, ho(a) — » 5(a) when a — > oo. We then have: £o = 1, £i = 0, 
£ 2 = 1/a, £3 = 2/a 2 , ^ 4 = ^ + Jj, £ 5 = 1| + If , etc. For a weakly polydisperse 
system we retain only the dominant terms of (3) in a 1/a expansion. From (3) we 
obtain then: Ai = ai(oo)£ 2 + 0(l/a 2 ), A 2 = ai(oo)£ 3 + a 2 (oo)(£ 4 - (g) + 0(l/a 3 ) = 
{ai(oo) + a2(oo)}^3 + 0(l/a 3 ), A3 = 01(00)^4 + 0(l/a 3 ), etc, where aj(oo) are the 
values of ai for a — > 00. Using the vdW expression (4) to evaluate (1) one finds, 
for ex. for t = 1.0 and tjq = 0.5, ai(oo) = 1.75 and 02(00) = —2.68. Using the 
corresponding numerical solutions found for h\(a) and /12(c) ( see Fig.l) it can be seen 
from Fig.2 that Ai/£ 2 ~ 1-75, A 2 /£ 3 « -0.93 and A 3 /£ 4 1.75 are obeyed to within 
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ten percent for a larger than, respectively, 40, 80 and 150. We can conclude thus 
that the weak polydispersity expansion (3) is valid (to dominant order) for Schulz 



These values do of course depend on the thermodynamic state but the case considered 
here (t = 1, 770 = 0.5) is representative of other t, 770 values. Note also that we have 
verified numerically that the radius of convergence of (3) with respect to 1/a is fairly 
sensitive to the total amount of polydispersity present. Allowing, for instance, the 
amplitude eo of the pair potential V(r; a, a') to depend on a and a 1 does reduce the 
radius of convergence of (3) considerably. From the above it follows that, ^ follows 
the universal law, ^ = ||, put forward in ref.4 whereas ^ follows the non- universal 
law, 4 s - = {1 + a2 l°°l || 3 -. We have verified that similar results can be obtained for 

' Ai L di (00) J £2 

different ho(a) distributions. Taking, for instance, a Gaussian for ho(a) similar results 
are found although £3 = for this case. This invalidates the conclusion of ref.4 that a 
particular importance should be attached to the skewness of ho(a). In conclusion, the 
general moment relation (3) can yield useful information about the phase behavior of 
weakly polydisperse systems but this information is in general not universal. 



distributions /10(c) with a dispersion 




smaller than, say, 0.1 (a « 100). 
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Figure Captions 



FIG. 1. The polydispersity distributions h n (a) of the parent phase (n = 0: full curve) 
(a Schulz distribution with the width parameter a = 50), the low-density (n=l: dotted 
curve) and the high-density (n=2: circles) daughter phases, as obtained by numerically 
solving the coexistence conditions of the van der Waals model of eq.(4) for t = 1, 
r/o = 0.5. The corresponding dimensionless densities of the coexisting daughter phases 
are r\\ = 0.106, 7/2 = 0.521 whereas for the monodisperse system one has, r/i = 0.103, 
r/2 = 0.608. Also shown are h\(a) — ho(a) (dashed curve) and [/12(c) — ho(a)] ■ 50 
(triangles) . 

FIG. 2. The ratio Afc/^ +1 (k = 1, 2, 3) versus \/a as obtained from the numerical so- 
lution of the van der Waals model of eq.(4) for t = 1, 770 = 0.5 and a Schulz distribution 
for /io(c). The symbols are as follows: circles(k=l), squares(k=2) and triangles (k=3). 
The dotted lines indicate their asymptotic (a — ► 00) values. The arrows indicate for 
each case the radius of convergence of the weak polydispersity expansion of eq.(3). 
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Fig.l, Xu et al., P.R.E. 
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Fig. 2, Xu et al., P.R.E. 
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